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CALCULATION  AND  MEASUREMENT  OF  SEPARATED  TURBULENT  BOUNDARY  LAYERS 

by 

P.  D.  Smith 
R.  C.  Hastings 
B.  R.  Williams 


SUMMARY 


An  inverse  integral  prediction  method  for  the  development  of  separated 
turbulent  boundary  layers  developed  from  the  lag-entrainment  method  is  des¬ 
cribed.  The  inverse  method  uses  the  concept  of  equilibrium  separated 
boundary  layer  flows  and  the  predicted  characteristics  of  such  flows  will  be 
compared  with  measurements  which  represent  the  first  known  demonstration 
that  equilibrium  separated  boundary  layers  can  be  realised  experimentally. 

In  these  experiments  the  data  were  obtained  with  a  single-component 
laser  Doppler  anemometer  usually  set  up  to  measure  streamwise  components  of 
mean  velocity  and  turbulence  intensity;  in  addition,  however,  one  pair  of 
profiles  of  the  mean  velocity  and  turbulence  intensity  normal  to  the  wall 
was  obtained.  Shear  stresses  were  not  measured  but  values  estimated  using 
the  data  in  conjunction  with  the  differential  momentum  equation  are  shown 
and  compared  with  the  patterns  already  revealed  by  earlier  work  on  attached 
equilibrium  flows. 

^  •'The  separated  flow  on  a  NACA  4412  aerofoil  has  been  measured  by 
Wadcock  using  the  flying  hot-wire  technique ._  The  measured  displacement 
thickness  is  used  as  input  to  the  inverse  prediction  method.  j£t  is  shown 
that  predicted  values  of  momentum  thickness  agree  with  the  measured  values 
but  that  the  calculation  predicts  a  pressure  rise  in  the  separated  region 
whereas  the  pressure  is  almost  constant  in  the  experiment.  The  result  of 
introducing  second  order  effects  into  the  calculation  is  shown.  The  equiva¬ 
lent  inviscid  flow  is  constructed  and  the  matching  of  the  equivalent  and 
real  flows  is  considered,. 
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!  THE  INVERSE  MODE  OF  THE  LAC-ENTRAINMENT  METHOD 

The  standard  direct  mode  of  the  lag-entrainment  method  computes  the  boundary  layer 
development  for  given  starting  conditions  and  a  specified  externally  imposed  pressure 
gradient.  East,  Smith  and  Merryman*  have  developed  an  inverse  mode  of  this  process  in 
which  the  pressure  distribution  is  predicted  for  a  given  displacement  surface.  To  do  this 
they  inverted  the  standard  equations  to  obtain  an  expression  for  the  pressure  gradient 
parameter  (0/Ue)  (dlWdx)  as  a  function  of  the  gradient  of  the  displacement  thickness, 
d-5*/dx  and  other  local  parameters,  ie 


6  dUe  =  J_  f  dSjl 

U  dx  F,  \  dx  1 

e  l 


a) 


where  F 


H  —  +  ( I  +0.2rM2) 


(s  -  “>  T4) 


dH 

dH, 


and 


-  H(H+2-M2)  +  (1  +0.2rM2)(H+  DH,  ~  +  0.4rM2(l  +0.2M2)(H+  1)  , 


6* 

T  ’  Hi 


(6-6*) 
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H  +  1  =  (H  +  1 )  ( 1  +  0.2M  )  , 


C£  is  the  skin  friction  coefficient, 

C£  is  the  entrainment  coefficient. 

As  written  by  East  et  al,  equation  (1)  is  singular  when  (dH,/dH)  -  0  but  this  is  easily 
overcome  by  writing  (1)  as 


l!!e  _  j/dH.  d^  _,\ 
ue  dx  "  ^2  v'dlT  dx  “  * \) 


(2) 


dH  dH. 

where  F’  =  — —  F„  and  Fj  =  -  F.  . 

1  dH  1  1  dH  1 

2 

Lock  has  pointed  out  that  for  interactive  calculations  involving  an  iteration 
between  a  boundary  layer  calculation  and  an  inviscid  flow  calculation  the  use  of  a  source 
strength  Z  rather  than  6*  is  to  be  preferred. 


Here 


,  d(pUe«*) 
p  dx 


te 


d6* 

dx 


so  that 
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Equilibrium  flows  are  defined  in  the  lag-entrainment  method  as  flows  in  which  the 
shape  of  the  velocity  and  shear  stress  profiles  in  the  boundary  layer  do  not  vary  with  x. 
The  parameters  used  to  characterise  compressible  equilibrium  flows  are 


C  U  dx 
f  e 


which  are  related  by  an  equilibrium  locus  of  the  form 


G 


6.432  (1  +  O.Bn) ( 1  +  0 


,04M2)] 


(4) 


Equation  (4)  may  be  expressed  in  a  form  more  suitable  for  separated  flows  (negative  C^  ) 

as 


0.02417  1  „  „  it 

- y  -  — ~  -  U  .0  —S' 

O*0.04M)  G  G 


(5) 


The  growth  of  the  equilibrium  boundary  layer  can  be  calculated  by  solving  the  momentum 
integral  equation  for  a  velocity  distribution  of  the  form 


U  “ 
e 


m 


x 


which  defines  the  pressure  gradient  parameter,  m  ,  as 


m 


U  dx  * 
e 


(6) 


(7) 


The  momentum  integral  equation  for  two-dimensional  compressible  flow  is 


de 

dx 


+ 


(H  +  2  -  M2)  i - 
e 


dU 

e 

dx 


hence  for  equilibrium  flow 


(H+  2  -M2) 


m8 

x 


1 


(8) 


(9) 


As  H  and  Cj  are  assumed  constant,  equation  (9)  can  be  solved  analytically  if  M 
is  also  assumed  independent  of  x  .  Strictly  such  an  assumption  is  not  consistent  with 
equation  (6)  but  it  is  a  valid  approximation  over  a  limited  length  of  boundary-layer 
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2  .  2 

development  during  which  the  change  in  M  is  small  compared  with  (H +  2  -  M  ).  Integrat¬ 
ing  equation  (9)  gives 


x  2^1  +  m(H  +  2  -M  )J 

Hence  for  these  flows  0  and  6*  tsince  H  is  constant)  vary  linearly  with  x  .  Sub¬ 
stituting  equations  (10)  and  (7)  into  equation  (5)  together  with  an  appropriate  skin  fric¬ 
tion  relationship  enables  an  explicit  expression  for  m  to  be  obtained  as  a  function  of 
R  and  H  .  This  function  is  portrayed  in  Fig  I .  In  general  for  any  given  value  of  m 
two  equilibrium  flows  are  possible,  one  separated  and  the  other  attached  although  for 
m  <  -0.2  two  attached  solutions  occur. 

3  MEASUREMENTS  IN  AN  EQUILIBRIUM  SEPARATED  BOUNDARY  LAYER 


In  order  that  predicted  and  measured  characteristics  of  equilibrium  separated  flow 

could  be  compared,  an  experiment  was  set  up  in  the  boundary  layer  tunnel  at  RAE  Bedford. 

In  this  tunnel  the  test  boundary  layer  is  the  one  which  grows  on  the  flat  ceiling  of  the 

working  section  and  it  can  be  subjected  to  a  variety  of  pressure  gradients  by  altering  the 

shape  of  the  flexible  floor  of  the  tunnel.  Equilibrium  attached  boundary  layers  had  pre- 

3 

viously  been  generated  and  studied  in  this  way  in  the  tunnel. 

In  this  instance,  the  requirement  was  for  a  conf iguration  that  would  permit  a 
separated  boundary  layer,  whose  external  velocity  (Ue)  varied  with  distance  according  to 
a  negative  power  law,  to  form  along  part  of  the  ceiling  at  the  same  time  as  there  was  a 
much  thinner  attached  boundary  layer  on  the  floor.  This  was  achieved  with  the  arrangement 
shown  diagrammatically  in  Fig  2.  Vortex  generators  were  used  to  inhibit  the  growth  of  the 
boundary  layers  on  the  tunnel's  side  walls  and  a  bleed  slot  was  inserted  at  a  joint  in  the 
tunnel  floor  to  thin  the  floor  boundary  layer. 

The  main  instrumentation  consisted  of  static  pressure  tappings  along  the  centreline 
of  the  ceiling  and  a  movable,  single-channel,  laser  Doppler  anemometer  to  measure  velocity 
components  parallel  to  the  ceiling.  The  anemometer  incorporated  Bragg  cells  to  introduce 
frequency  shifts  into  its  light  beams  so  that  the  signs  of  the  measured  velocity  compo¬ 
nents  could  be  determined.  In  addition,  the  anemometer  was  temporarily  adapted  to  measure 
vertical  velocity  components  at  one  streamwise  station  along  the  working  section. 

The  streamwise  variation  of  pressure  measured  along  the  centreline  of  the  ceiling 
is  shown  in  Fig  3a  and  the  extent  of  the  ceiling  and  floor  boundary  layers  in  Fig  3b. 

The  region  of  equilibrium  flow,  between  an  upstream  zone  of  development  toward  equilibrium 
and  a  downstream  zone  in  which  the  ceiling  and  floor  boundary  layers  are  believed  to  have 
been  sufficiently  close  to  interfere  with  each  other  is  shown  in  Fig  3a. 

In  the  equilibrium  region,  the  pressure  variation  corresponds  to 


U  «  x 
e 


where  m  »  -0.122  if  x  is  measured  from  the  virtual  origin  of  the  equilibrium 

boundary  layer. 
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The  displacement  and  momentum  thicknesses  (<5*,0)  for  the  test  boundary  layer  on  the 
centreline  of  the  ceiling  are  shown  in  Fig  4.  This  shows  that  the  linear  growth  charac¬ 
teristic  of  equilibrium  flow  was  achieved.  Measurements  offset  200  mm  on  either  side  of 
the  centreline  were  also  made.  The  offset  momentum  and  displacement  thicknesses  agree 
with  centreline  values  within  about  8%  and  10%  respectively.  The  flow  appears  to  be 
slightly  asymmetrical  with  values  to  starboard  agreeing  better  with  the  centreline  values 
than  those  to  port;  the  latter  correspond  to  a  slightly  thinner  boundary  layer  whose 
shape  parameter  H  is  between  1%  and  10%  lower  than  comparable  centreline  values. 

The  centreline  measurements  are  consistent  with  the  following  boundary  layer  growth 
properties: 

=  0.0193  ,  =  0.104  • 

dx  ’  dx 

Consequently  mean  values  of  H  =  6*/0  and  J  =  (H-l)/H  are  5.38  and  0.814  respectively. 

The  experimental  values  J  =  0.814  and  m  =  -0.122  may  be  compared  with  the  pre¬ 
dicted  (m,J)  locus  discussed  in  section  2,  for  which  m  =  -0.122  gives  J  =  0.839  .  It 
will  be  seen  that  J  does  not  fall  on  but  below  the  predicted  locus;  in  fact  it  is 
roughly  as  close  to  the  locus  as  were  the  attached  flow  measurements  of  Ref  3. 

Centreline  profiles  of  streamwise  velocity  ratio  U/Ue  at  a  series  of  streamwise 
locations  are  shown  in  Fig  5.  They  exhibit  an  encouraging  degree  of  similarity,  though 
it  has  to  be  admitted  that  the  streamwise  length  covered  by  the  measurements  is  rather 
small . 

Profiles  near  the  wall  are  shown  to  expanded  scales  in  Fig  6.  There  is  no  evidence 
that  the  universal  law  of  the  wall  for  attached  flow,  with  an  appropriately  modified 
definition  (-tw/o)^  of  ,  can  be  applied  to  the  reverse  flow  profile.  This  is 
because  the  maximum  reversed  velocities  occur  very  close  to  the  surface. 

Fig  7  shows  the  only  profile  of  the  vertical  velocity  ratio  V/Ug  that  was 
measured.  The  profile  shape  is  generally  consistent  with  two-dimensional  equilibrium 
flow  if  U/Ue  behaves  as  shown  in  Fig  5.  However  values  of  V/Ug  are  rather  too  large 
in  general. 

~2  2 

Normalised,  mean  square,  vertical  and  horizontal  components  of  turbulence,  u  /U 

T  2  e 

and  v  /Ug  ,  are  shown  in  the  next  figure  (Fig  8)  for  the  one  station  where  both  are 

available.  The  distributions  of  these  quantities  may  be  compared  with  the  corresponding 

distributions  for  a  flow  near  separation  (Flow  7)  from  Ref  3.  Both  components  tend  to 

peak  near  y/6  »  0.6  for  both  flows,  but  the  peaks  for  both  components  are  rather  more 

than  twice  as  high  in  the  separated  flow  as  in  the  attached  flow. 

The  ratio 


for  the  separa*  ed  flow. 
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This  indicates  that  normal  Reynolds  stress  effects  are  considerably  more  important 
relative  to  shear  stress  effects  in  the  momentum  integral  equation  than  would  be  expected 
on  the  basis  of  results  for  equilibrium  attached  flow.  It  is  found  for  the  latter  in 
Ref  3  that  the  ratio  is  always  about  0.072 J  so  that  the  expected  value  would  be  about 
0.059  for  J  =  0.814  . 

Since  the  shear  stress  distribution  across  the  separated  boundary  layer  was  not 
measured,  it  was  thought  worthwhile  to  estimate  it.  This  was  done  in  the  following  way. 

It  was  assumed  that  viscous  effects  would  be  important  only  in  a  tiny  region  close  to  the 
surface  and  could  be  neglected.  On  this  premise  it  was  further  assumed  that:  the  velocity 
profiles  at  different  streamwise  locations  would  have  the  same  shape  (that  is  U/U^  would 
be  a  function  of  y/6  only),  6  a  x  ,  and  a  x  ,  if  x  is  measured  from  the  virtual 
origin  of  the  equilibrium  layer.  These  assumptions  enable  the  differential  equation  of 
boundary  layc  motion  to  be  reduced  to  an  ordinary  differential  equation  if  it  is  also 
assumed  that  the  shear  stress  term  -UV  is  expressible  as  the  product  of  a  function  of 
y/6  and  a  factor  proportional  to  x^m  .  The  differential  equation  developed  in  this  way 
was  integrated  numerically  to  find  the  shear  distribution,  with  appropriate  experimental 
information  fed  in,  and  the  starting  condition  that  the  shear  was  zero  at  the  surface. 

This  condition  had  the  advantage  that  no  assumed  value  of  skin  friction  coefficient  (C.) 

2  1 

was  required  to  start  the  calculation  but  meant  that  the  computed  estimate  of  x/pU^ 
should  tend  to  -Cj/2  rather  than  zero  at  the  edge  of  the  boundary  layer. 

Two  basic  calculations  were  performed:  a  first  approximation  to  the  shear  distribu¬ 
tion  was  obtained  Ly  assuming  that  static  pressure  was  constant  across  the  thickness  of 
the  boundary  layer,  and  a  second  approximation  in  which  3P/3y  =  -(3pV  /3y)  and  the 
effect  of  the  term  3U^/3x  on  the  streamwise  momentum  balance  was  included.  Results  are 
shown  in  Fig  9.  The  first  approximation  yielded  a  shear  stress  distribution  with  a  rela¬ 
tively  high  peak  and  thereafter  a  fairly  sharp  decay.  The  value  outside  the  boundary 
layer  was  negative  and  not  the  small  positive  value  to  be  expected.  Taking  a  second 
approximation  produced  a  slightly  lower  peak  followed  by  a  more  gradual  and  plausible 
decay;  however  values  outside  the  boundary  layer  appeared  to  be  too  high. 

An  additional  case  was  computed  in  which  the  input  values  of  36/3x  and  |m|  were 
respectively  increased  and  decreased  by  2.7%.  This  case,  too,  is  shown  in  Fig  9.  The 
residual  shear  at  the  outer  edge  of  the  boundary  layer  was  halved.  The  2.7%  changes 
illustrated  had  no  appreciable  effect  on  the  shear  stress  estimate  for  the  inner  two- 
thirds  of  the  boundary  layer.  The  estimate  for  the  inner  third  is  almost  linear  and 
clearly  the  relation  3P/3x  =  3x/3y  ,  valid  at  y  =  0  ,  is  the  dominant  influence  on  the 
profile  in  that  region. 

The  next  figure,  Fig  10,  shows  some  shear  stress  distributions  taken  from  the 
attached-flow  measurements  of  Ref  3.  The  figure  shows  how  the  maximum  normalised  shear 
stress  both  moves  away  from  the  surface  and  increases  in  magnitude  as  separation  is 
approached.  The  shear  stress  estimates  for  the  separated  flow  are  plausible  in  general 
terms,  in  that  they  indicate  that  the  trend  continues. 
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All  in  all,  the  experiment  on  the  one  separated  flow  with  an  adverse  pressure 
gradient,  which  it  has  so  far  been  possible  to  set  up,  confirms  that  the  concept  of  equi¬ 
librium  separated  flow  is  valid;  however  detailed  examination  of  the  measurements  indicate 
that  the  flow  studied  is  not  truly  in  two-dimensional  equilibrium.  The  profile  of  V/Ug 
is  not  completely  consistent  with  the  U/Ug  profile  and  the  gradient  of  the  latter,  out¬ 
side  the  boundary  layer  (this  may  just  be  discerned  in  Fig  5)  is  rather  too  large  to  be 
compatible  with  the  growth  of  the  displacement  thickness.  Furthermore,  the  velocity  com¬ 
ponent  profiles,  surface  pressure  gradient  and  boundary  layer  growth  do  not  lead  to  com¬ 
pletely  satisfactory  estimates  of  shear  stress.  The  importance  of  second-order  effects  in 
this  boundary  layer  is  noteworthy;  it  may  be  that  third  order  effects  are  not  completely 
negligible. 

4  CALCULATIONS  OF  SEPARATED  FLOW  ON  A  TWO-DIMENSIONAL  AEROFOIL 

The  use  of  the  inverse  boundary-layer  method  developed  by  East  et  al}  is  demonstrated 
by  calculating  the  development  of  the  separated  flow  at  the  trailing  edge  of  a  NACA  4412 

4 

aerofoil,  using  the  measurements  of  Wadcock  for  the  boundary  conditions.  Having  demon¬ 
strated  the  ability  of  the  inverse  boundary  layer  method  to  calculate  this  class  of  flow, 
the  flow  about  the  complete  aerofoil  is  calculated  by  matching  an  external  inviscid  incom¬ 
pressible  flow  with  the  inverse  boundary  layer  by  the  scheme  developed  by  Le  Balleur"’  and 

2 

further  described  by  Lock  .  It  is  shown  that  the  solutions  converge  for  examples  with 
separating  and  separated  flow  at  the  trailing  edge. 

4. 1  Calculation  of  the  development  of  the  boundary  layer 

Measurements  by  a  flying  hot-wire  of  the  separated  flow  on  a  NACA  4412  aerofoil  at 
o  4 

13.87  incidence  have  been  reported  by  Wadcock  .  The  displacement  thickness  measured  by 
Wadcock  has  been  used  as  input  to  the  inverse  lag  entrainment  method.  It  is  shown  that 
the  predicted  values  of  momentum  thickness  and  shape  factor  agree  with  the  measured 
values,  but  the  predicted  values  of  the  velocity  do  not  agree  with  the  measured  values. 

It  is  demonstrated  that  quite  large  changes  in  the  displacement  thickness  produce  only 
small  changes  in  the  predicted  values,  thus  the  disagreement  cannot  be  explained  by 
errors  in  the  measurements  in  the  experiment.  However  a  change  in  the  H  ~  Hj  relation¬ 
ship  produces  a  marked  change  in  the  predicted  velocities,  which  suggests  that  some 
future  work  is  required  to  derive  an  appropriate  H  ~  Hj  relationship.  To  allow  an 

exploration  of  a  fuller  range  of  H  ~  H.  relationship  the  form  of  the  lag  entrainment 

2 

equations  suggested  by  Lock  is  used  so  they  can  be  integrated  through  the  singularity 
which  occurs  in  the  form  given  by  East  et  al *  if  dHj/dH  **  0  .  It  should  be  noted  that 
the  predicted  velocities  must  be  corrected  for  second  order  effects  before  they  can  be 
compared  with  the  measured  velocities.  East^  has  pointed  out  that  in  a  proper  second 
order  treatment  of  the  momentum  integral  equation  the  velocity  distribution  from  the 
equivalent  inviscid  flow  should  be  used  rather  than  the  distribution  from  the  real  viscous 
flow.  By  the  converse  argument  it  follows  that  the  inverse  boundary  layer  equation  will 
predict  the  velocity  in  the  equivalent  inviscid  flow.  At  separation  the  second  order 
terms  are  important  which  will  give  a  significant  difference  between  the  equivalent  and 
real  flows.  An  attempt  is  made  to  quantify  these  differences  and  perform  a  valid  compari¬ 
son  of  the  measured  and  predicted  velocity  distributions. 
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The  experimental  data  used  in  these  calculations  are  given  in  Table  1 .  The  surface 
velocities  are  interpolated  from  the  data  in  Wadcock's  PhD  thesis  and  the  source  strength 
is  calculated  by  numerical  differentiation.  In  Fig  1 1  the  development  of  H  is  shown  to 
compare  well  with  the  experimental  values  and  an  equally  good  agreement  was  obviously 
found  for  the  momentum  thickness.  These  calculations  were  performed  with  the  H  ~  Hj 
relationship  given  by  East  et  al)  , 


H  =  4.5455  +  295  exp(-  3.325H) 

dH  exp(3.325H) 

dH,  =  "  980 


(11) 


-  2 

Results  are  also  shown  for  a  ’modified’  East  H  -  H|  relationship  in  which  0.01  (H  -  1 ) 
has  been  subtracted  from  the  expression  for  H,  .  This  results  in  a  reduction  in  the 
coefficient,  dH/dH(  in  the  differential  equation  for  H  and  consequentially  a  slower 
growth  of  H  . 

The  velocity  distribution  predicted  by  the  inverse  boundary  layer  method  is  given 
in  Fig  12  for  both  H  -  H,  relationships.  There  is  a  marked  disagreement  between  the 
predicted  and  measured  values,  with  the  predicted  velocities  decreasing  throughout  the 
separated  region. 

It  is  possible  that  the  displacement  thickness  was  measured  incorrectly.  However 
an  increase  in  the  displacement  thickness  by  35%,  which  should  exceed  any  experimental 
error,  only  increased  the  velocity  at  the  trailing  edge  by  0.02  whilst  the  difference 
between  equivalent  and  real  flows  is  0.135  as  shown  in  Fig  12. 

2 

By  integrating  the  normal  momentum  equation  Lock  shows  that,  to  second  order,  the 
difference  between  the  pressure  at  the  wall  in  the  equivalent  inviscid  flow,  P^  and 
the  pressure  in  the  real  viscous  flow,  P^  is  given  by 


P. 

1W 


-  P 

w 


p ,  U?  K*(8  +  6) 
lw  1W 


( I  2) 


where  K*  ij  the  curvature  of  the  displacement  surface.  In  Fig  13  this  correction  has 
been  added  to  the  equivalent  inviscid  flow  and,  bearing  in  mind  the  difficulty  of  evalua¬ 
ting  the  curvature  accurately,  there  is  reasonably  good  agreement  with  the  real  viscous 
flow  up  to  separation.  After  separation  the  displacement  surface  is  not  curved  and  this 
correction  has  a  negligible  effect. 

The  effect  of  using  different  H  -  H.  relationships  was  examined  and  the  relation- 

.  .  .  '  7 

ships  are  plotted  in  Fig  14  and  are  listed  below:  Green 


1 .09289 

+  n  s 

p-  n 

1.09289 

Vh  -  \) 

*»■  U  .3 

I1.12J 

> 


-  -  1  .85549(H- 1)~2,09289  +  0.48279(H- l)0,09289  .  J 


(13) 
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Hast 


1 


=  3.15  + 


dH 

dH, 


(H-  1) 
I  .72 


1.72 
(H-  ! ) 

2 


H  «  1  .6 


(14) 


As  equation  (11)  if  H  >  1 .6  . 

The  equations  for  the  inverse  method  for  incompressible  flow  take  the  form: 


g  -  -£-«•«*§ 


£  du 

U  dx 


(H 


1  Cf 
+  1)  2 


C'  CE  -  H.£ 


dH 

dx 


1  <V  -  HV 


dC 


E 

dx 


*  F(v  \) 


2.8 


H  +  H j  \(C t)EQq  "  >C\ 


where  the  last  equation  is  the  'lag  equation'  for  C£  , 

dH, 


D  =  Hi  -  H  IF 


/e  du\  e_du 

\V  dx/EQ  ‘  U  dx 


(15) 


(16) 


(17) 


(18) 


which  is  never  negative,  and  now  the  source  strength  Z  is  regarded  as  'known'  and 
(1 /U) (dU/dx)  as  'unknown'.  It  will  be  seen  in  section  5  that  the  source  strength  rather 
than  the  displacement  thickness  ties  in  with  the  scheme  for  matching  the  viscous  and 
inviscid  flows. 

The  results  for  the  different  H  —  H  relationships  are  indicated  in  Figs  11  and  12. 
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The  expression  originally  used  by  Green  produces  a  marked  change  in  both  the  velocity  and 
the  development  of  H  ;  but  neither  of  these  changes  appears  very  realistic  and  more 
experimental  evidence  is  required  to  determine  an  appropriate  H  ~  H|  relationship. 

4 . 2  Matching  of  the  inviscid  and  viscous  flows 

In  the  previous  section  it  has  been  demonstrated  that  the  boundary  layer  equations 
can  be  integrated  through  the  singularity  at  separation  if  the  equations  are  cast  in 
inverse  form  for  which  prescribed  values  of  the  source  strength  Z  predict  the  velocity 
gradient  parameter  (1/U) (dU/dx)  (see  equation  (16)).  It  seems  natural  to  use  the  solu¬ 
tion  of  the  inverse  boundary  layer  equations  as  input  to  an  inverse  inviscid  method  to 
obtain  a  new  estimate  of  the  displacement  surface  ( ie  Z) .  However  it  may  be  difficult  to 

g 

solve  the  inviscid  problem  since  Motz  has  indicated  that  there  is  a  singularity  in  the 
inviscid  flow  at  the  point  where  the  boundary  conditions  change  from  direct  (shape 
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specified)  to  inverse  (velocity  specified),  ie.  at  the  transition  from  calculations  with 

the  direct  boundary  layer  method  to  the  inverse  form.  This  difficulty  is  overcome  if  the 
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'semi-inverse'  technique  as  described  by  Carter  ,  Le  Balleur  ,  Lock',  Wigton  and  Holt 
is  adopted. 

Such  a  technique  is  illustrated  in  Fig  15.  An  estimate  of  the  source  strength  on 
the  aerofoil  and  the  sink  strength  in  the  wake  is  used  to  calculate  (a)  the  inviscid 
velocity  U.  by  a  direct  method,  and  (b)  a  velocity  U.  from  the  inverse  boundary 

iw  IW 

layer  method.  The  two  velocities  are  compared  and  if  the  scheme  has  not  converged  then 
the  source  strengths  are  updated  by  the  formula 


.n+1 


V 

Uiw 


dUV 

lVi 

dx 


6  dU^ 

6  iw 


-|(n) 


U1 2 * 4 

iw 


dx 
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The  analysis  of  Le  Balleur  treats  the  error  I  -  I  as  a  superposition  of  Fourier 
modes  and  this  leads  to  an  expression  for  the  relaxation  factor  w  ,  which  for  incom¬ 
pressible  flow  reduces  to 

B 

W  ~  1  +  Ben/ Ax 


where  Ax  is  the  local  step  length  and  B  =  (H  +  I)(Hj  -  HdH/dH( )  / (dH/dtlj )  which  is  the 
coefficient  of  (0/U)(dU/dx)  in  the  expression  for  the  source  strength  derived  from 
equation  (16). 

The  direct  inviscid  solution  for  incompressible  flow  can  be  obtained  by  a  panel 

method  and  in  this  case  the  symmetric-singularity  model  developed  by  Newling*'  is  used. 

The  initial  development  of  the  boundary  layer  is  laminar  and  this  is  calculated  directly 

12.. 

by  the  method  due  to  Thwaites  .  The  transition  to  a  turbulent  boundary  layer  can  occur 

1  3 

through  natural  transition  which  is  predicted  by  Granvilles'  criterion  ,  or  a  laminar 

1 4 

separation  bubble  predicted  by  Horton's  method  or  forced  by  a  transition  trip.  The 
development  of  the  turbulent  boundary  layer  is  calculated  in  either  direct  or  inverse 
mode  by  the  lag  entrainment  method.  The  turbulent  boundary  layer  calculation  can  be 
switched  from  direct  to  inverse  mode  by  any  pre-determined  criterion:  for  the  following 
calculations  the  switch  occurs  when  H  =  1.7  .  The  development  of  the  combined  wakes 
from  the  upper  and  lower  surfaces  is  also  calculates  by  the  lag  entrainment  method,  using 
the  inverse  mode.  In  the  present  calculations  second  order  effects  have  not  been 

included  in  either  the  calculation  of  the  boundary  layer  or  the  construction  of  the 

2 

equivalent  inviscid  flow.  However  following  the  method  described  by  Lock  these  effects 
can  be  included. 


The  method  has  been  used  to  calculate  the  flow  about  a  NACA  4412  aerofoil  at  a 
Reynolds  number  of  1.5  *  10^.  Transition  was  forced  on  the  upper  and  lower  surfaces  at 
x/c  «  0.025  and  x/c  «  0.103  respectively  so  that  the  results  could  be  compared  with  the 

4 

experiment  described  by  Wadcock  . 
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The  first  calculation  was  performed  at  an  angle  of  incidence  of  8°  for  which  the 
flow  is  almost  separating  at  the  trailing  edge.  In  Fig  16  it  can  be  seen  that  the  scheme 
converges  to  a  lift  coefficient  slightly  above  the  experimental  value  for  a  calculation 
with  the  semi-inverse  technique  applied  on  the  upper  surface  wherever  H  >  1.7  and  in 
the  wake.  A  convergent  solution  is  also  produced  if  the  semi-inverse  technique  is  used 
for  all  the  region  occupied  by  the  turbulent  boundary  layer  on  the  upper  surface.  In 
Fig  17  the  pressure  distributions  are  compared  and  the  theoretical  results  have  a  steeper 
adverse  pressure  gradient  on  the  suction  surface.  In  the  previous  section  it  was  shown 
that  a  correction  of  the  equivalent  inviscid  flow  for  curvature  effects  would  reduce  this 
difference.  It  should  also  be  noted  that  passive  flow  control  was  added  to  Wadcock's 
experiment  to  obtain  two-dimensional  flow  at  14°  incidence.  The  same  control  was  used  at 
8°  incidence  so  the  flow  may  not  have  been  completely  two-dimensional. 

At  an  incidence  of  14°  the  method  again  converges  smoothly  (see  Fig  18)  but  the  lift 
is  now  5%  above  the  experimental  value.  For  the  first  nine  iterations  the  flow  is 
attached  to  the  aerofoil  but  in  Fig  19  it  can  be  seen  that  the  separation  point  moves 
forward  as  the  number  of  iterations  increases  until  it  reaches  x/c  =  0.918.  From  Fig  20 
the  separation  point  in  the  experiment  is  further  forward  and  most  likely  lies  between 
x/c  =  0.70  and  x/c  =  0.75.  The  theoretical  method  again  predicts  an  adverse  pressure 
gradient  in  the  separated  region  and  the  difference  in  the  pressure  distribution  in  the 
separated  region  are  a  major  contributor  to  the  discrepancy  in  shown  in  Fig  18. 

The  pressure  distribution  in  the  leading-edge  region  for  the  theory  and  the  experi¬ 
ment  are  plotted  in  more  detail  in  Fig  21 .  In  the  theoretical  calculation  there  is  a 

laminar  separation  bubble  which  reattaches  ahead  of  the  position  of  the  transition  trip 

4  ... 

in  the  experiment.  Wadcock  did  not  detect  any  laminar  separation  bubble  in  the  experi¬ 
ment.  At  the  transition  trip  no  increase  in  momentum  thickness  has  been  included  in  the 
calculation  and  this  may  explain  the  discrepancy  in  the  development  of  momentum  thickness 
shown  in  Fig  22.  In  Fig  23  there  is  a  larger  discrepancy  in  H  and  this  behaviour  ties 
in  with  the  observed  and  predicted  separation  positions. 

The  semi-inverse  scheme  provides  a  method  of  obtaining  converged  solutions  for  flows 
with  small  regions  of  separated  flow  at  the  trailing  edge.  However  the  details  of  the 
calculation  have  not  been  in  close  agreement  with  the  experiment  of  Wadcock.  Although 
there  is  a  reasonable  agreement  between  the  pressure  distributions  over  the  first  70%  of 
the  aerofoil  there  is  not  such  good  agreement  in  the  boundary  layer  parameters  over  the 
same  region.  In  the  separated  region  there  is  a  discrepancy  in  the  measured  and  predicted 
pressures  which  was  discussed  in  section  4.1. 

In  the  future  second  order  effects  will  be  included  in  the  boundary  layer  calcula¬ 
tion  and  the  matching  of  the  inviscid  and  viscous  flows  and  modifications  to  the  H  —  Hj 
relationship  for  the  separated  flow  region  will  be  investigated. 
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Measured  and  predicted  velocity  distributions  for  NACA  4412 
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■ection  added  to  equivalent  inviscid  flow 
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Fig  15  Flow  diagram  for  semi-inverse  scheme 
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